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ABSTRACT 

A family of potential-density pairs has been found for spherical halos and bulges of 
galaxies in the Newtonian limit of scalar-tensor theories of gravity. The scalar field is 
described by a Klein-Gordon equation with a source that is coupled to the standard 
Poisson equation of Newtonian gravity. The net gravitational force is given by two 
contributions: the standard Newtonian potential plus a term stemming from massive 
scalar fields. General solutions have been found for spherical systems. In particular, 
we compute potential-density pairs of spherical galactic systems, and some other as- 
trophysical quantities that are relevant to generating initial conditions for spherical 
galaxy simulations. 

Key words: stellar dynamics - elementary particles - Galaxy: kinematics and dy- 
namics - Galaxy: halo- Galaxy: bulge - Galaxy: structure 



1 INTRODUCTION 

In recent years there has been important progress in under- 
standing the dynamics that led to the formation of galaxies. 
Two and three dimensional iV-body simulations of galax- 
ies and protogalaxies have been computed using millions 
of particles, giving a more realistic view of how galaxies, 
quasars, and black holes could have formed JBarnes 1998 
|Barnes fc Hernquist 1992| l. 

The composition of Universe at the time galaxies formed 
could theoretically be very varied, including visible and 
dark baryonic matter, non-baryonic dark matter, neutri- 
nos, and many cosmological relics stemming from symme- 
try breaking processes predicted by high energy physics 
(Kolb & Turner 1990J) . All these particles and fields, if 
present, should have played a role in structure forma- 
tion. Accordingly, recent independent observational data; 
CMBR at various angular scales (Ide B crnardis ~ al 2000| 
IBennett et al 200311 . type la supernovae (IRiess et al 19981 
IPerlmutter et al 19991 IRiess et al 20011 . as well as the 2dF 
Galaxy Redshift survey (IPeacock 20021 lEfstathiou 200211 . 
suggest that Q — Qa + O m w 1, or Qa ~ 0.7 and 
fl m ~ 0.3, implying the existence of dark energy and 
dark matter, respectively. One particular candidate for 
dark energy is a scalar field usually called quintessence 
jCaldwell, Dave fc Steinhard 1998| l . In this way galaxies 
are expected to possess dark components and, in accor- 
dance with the rotation curves of stars and gas around 
the centres of spirals, they might be in the form of 



halos (IQstriker fc P eebles 1973$ and must contribute to 
at least 3 to 10 times the mass of the visible matter 
HKolb fc Turner 19901 . 

In order to construct numerical galaxies in a consistent 
way using gravity and particle physics components, one 
must find the correct physical scheme. The theoretical 
framework to explain the existence of dark components 
finds its origin in theories of elementary particles physics, 
with the addition of the action of gravity. There exist many 
theories (grand unification schemes, string theories, brane 
worlds, etc) that involve such physics, but scalar-tensor 
theories (STT) of gravity are typically found to represent 
classical effective descriptions of such original theories 
l |Green, Schwarz fc Witten 1988^ . In this way, the scalar 
fields of these theories are the natural candidates to be 
the quintessence field l|Qfidw^h^Dave^^^tehiharci^)98| 
B oiss eau, Esposito-Farese^oku"skiT^St£TObTisky^000] 
kmendola 20011 . as a remnant of some cosmological func- 
tion that contributes Qa ~ 0.7 today. It has been even 
suggested that the quintessence field is the scalar field 
that also acts on local planetary scales ( |Fujii 20 00 ) or on 
galactic scales (Matos & Guzman 20011. Moreover, massive 
scalar fields might account to the dark matter components 
of galaxies in the form of halos. 

Motivated by the above arguments, we study in the 
present report some STT effects in galactic systems. Specifi- 
cally, we compute potential-density pairs coming from such 
theories in their Newtonian approximation. The results pre- 
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sented here will be useful in constructing numerical galaxies 
that are consistent with grand unification schemes. 

This paper is organised as follows: in section 2 we 
present the Newtonian approximation of a general scalar- 
tensor theory of gravity. Solutions are presented in terms of 
integrals of Green functions, and some expressions for the 
velocities and dispersions of stars in galaxies are given. In 
section 3, solutions for a family of potential-density pairs 
UDehnen 19931 and the Navarro, Frenk, and White (NFW) 
model INavarro et al 1996-971 are presented, and some rele- 
vant observational quantities are computed. In section 4 we 
present conclutions. 



V 2 ip = 47rp . (6) 
The new Newtonian potential is now given by 

$jv = -hoo =tp- i<£- (7) 

The next step is to find solutions for this new Newto- 
nian potential given a density profile, that is, to find the 
so-called potential-density pairs. First we present the solu- 
tions for point -like masses and second for general density 
distributions. 



2 SCALAR FIELDS AND THE NEWTONIAN 
APPROXIMATION 

Let us consider a typical scalar-tensor theory given by the 
following Lagrangian 



C = 



16tt 



-SR+^{dS) 2 -V{6) 



(l) 



from which we get the gravitational equations, 



1 LO 

87rT M „ + -Vg^v + — d^SdvS 
Z q> 



g^u us 



and the scalar field equation 
6V' -2V 1 



US , 
V 3 + 2a; 

where ()' = 



3 + 2oj 1 
4f,- In 



8ttT - J(dS) 2 



(2) 



(3) 



the past we have studied 



( Cervantes- Cota & Dehnen 1995), among others, this type 
of non-minimal coupling between gravity and matter fields 
in a cosmological context in the very early Universe, where 
all relativistic effects have to be taken into account. In the 
present study, however, we want to consider the influence 
of scalar fields within galaxies, and therefore we need to de- 
scribe the theory in its Newtonian approximation, that is, 
where gravity and the scalar fields are weak and velocities 
of stars are non-relativistic. We expect to have small devi- 
ations of the scalar field around the background field, de- 
fined here as (S) — G^ 1 = 1. If one defines the perturbation 
= S— 1, the Newtonian approximation gives ( Hclbig 1991 ) 

I, 

2 



Roo = -V 2 fooo 



4vrp _ I V 2 <; 



V S — m S = — Snap 

where we have introduced 
SV - 2V 



(4) 
(5) 



2 7 

= m <p • 



m 2 kS 2 + . 



3 + 2a; 

and a = 1/(3 + 2a;). In the above expansion we have set the 
cosmological constant term equal to zero, since on galactic 
scales its influence should be negligible. We only consider 
the influence of dark matter due to the boson field of mass 
m governed by Eq. 0, that is a Klein-Gordon equation 
with a source. 

Note that Eq. Q can be cast as a Poisson equation for 
^=(1/2)0100 + 0), 



2.1 Point— like masses 

The solution of these equations is known and given by 



s = 


2au x , 


hoo = 


— 2u — 2au\ 


where 




u — 


\ - m s 

^ |r-r.| ' 

S 


U\ = 


\ - m s 

^ |r-r s | 6 



(8) 



(9) 



with m s a source mass, and the total gravitational force on 
a particle of mass rrii is 



-V/loo = rata, 



(10) 



where A = h/mc is the Compton wavelength of the effec- 
tive mass (m) of some elementary particle (boson) given 
through the potential V(S) and u>(S). In what follows we 
will use A instead of m" 1 . This mass can have a range of 
values depending on particular particle physics models. The 
potential u is the Newtonian part and u\ is the dark mat- 
ter contribution which is of the Yukawa type. There are two 
limits: On the one hand, for if r 2> A (or A — > 0) one re- 
covers the Newtonian theory of gravity. On the other hand, 
for if r <C A (or A — > oo) one again obtains the Newto- 
nian theory, but now with a rescaled Newtonian constant, 
G — > Gjv(l + a). There are stringent constraints on the pos- 
sible A— a values determined by measurements on local scales 
HFischbach fc Talmadge 1999) . 

In the past the above solutions have been used to 
solve the missing mass problem in spirals USanders 19841 
Eckh ardTT993t as an alternative to considering a distribu- 
tion of dark matter. This was done assuming that most of 
the galactic mass is located in the galactic centre, and then 
considering the centre to be a point source. In our present 
investigation we do not avoid dark matter, since our model 
predicts that bosonic dark matter produces, through a scalar 
field associated to it, a modification of Newtonian gravity 
theory. This dark matter is presumably clumped in the form 
of dark halos. Therefore we will consider in what follows 
that a dark halo is spherically distributed along an observ- 
able spiral and beyond, having some density profile (section 
3). Next, we compute the potentials, and some astrophysical 
quantities, for general halo density distributions. 
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2.2 General density distributions 

General solutions to Eqs. JjsJ and (JSJ can be found in terms 
of the corresponding Green functions 



drsJ £^L + B.C., 



= 2n I dr / (rfl , )e ' r r3 ' A +B.C. 



|r - r s 

and the new Newtonian potential is 



(11) 
(12) 



dr. 



„ / „rj>^ ' ' A . ».c. 



(13) 



The first term of Eq. 1131 . given by ij), is the contribution 
of the usual Newtonian gravitation (without scalar fields), 
while information about the scalar field is contained in the 
second term, that is, arising from the influence function de- 
termined by the Klein-Gordon Green function, where the 
coupling Lu (a) enters as part of a source factor. 

Given that we are interested in spherical galaxies and 
bulges in what follows we will only consider the case of 
spherical symmetry. Additionally, we use flatness boundary 
conditions (B.C.) at infinity, such that the boundary terms 
in Eqs. 1111 - 112(1 are zero. Moreover, regularity conditions 
must be applied to spatial points where the potentials are 
singular. For spherical systems these conditions mean that 
dip/dr — d(j>/dr — at the origin. Accordingly, performing 
the integrals in these equations, we obtain ijj for a system of 
radius _R, 



and for 




^f 1 + o(r a ) ;r a <R a 
; r a > R a 



(14) 




X V„)+ Sinh(r r a/Xa) q(r a ) ;r a <R a 



(16) 



where we have defined the functions 

r- 2 

m(r a ) = / dx x p(x) , 
Jo 

o(r a ) = j dx xp{x) , 

. . _ f Ta , 2 sinh(:r/A a ) . 

p(r a ) = / dxx p(x) , 

Jo 

i \ _ f Ra j 2 exp(-x-/A a ) . 
q(r a ) = dxx -2 '-p(x) , 

J ra X / X * 
and we have written the density as 

Here M is a reference or scaling mass and a is a scaling 
length that we use to define the following quantities: r a = 
r/a, R a = R/a, \ a = A/a, and r] = a/X a . Then, the new 
Newtonian potential can be written as 



(17) 



$jv(r a ) = F(r„), 

where 



(18) 



F(Ta) 



. ainh(yA o) g(ra) 



r a /X a 



P{r a ) 



;r a < R a 



;r a > R a 



(19) 



2.3 Intrinsic and projected properties 

We now want to present the intrinsic and observable quanti- 
ties of interest that characterize the steady state of spherical 
galaxies and bulges. These include the circular velocity, ve- 
locity dispersion, distribution function, projected velocity 
dispersion, and projected density. 

The circular velocity for a spherical galactic system is, 
from the above equations, 



d<I>jv 
dr 



M dF(r a ) 

Ta ~a 

a dr a 



(20) 



For r a < R a one obtains straightforwardly the circular ve- 
locities for stars inside a dark halo of size R. For r a > -Ra 
one gets 



m(R a ) 



M 



where 



Q(r a ) 



a [o(R a ) sinh(i? a /A Q ) - Q(R a )] 



-(1 + ^): 



A„ 



ra/Aa A a ' 

dx x o(x) cosh(x/A a ) , 



(21) 



(22) 



which gives the circular velocities outside the halo (if there 
are stars there). From Eqs. 1191 and 1201 one obtains, after 
redefinition of the multiplicative constant, the expression 
found in the past for rotational velocities due to point masses 
ISanders 19841 . 



(15) 



A full description of stars in a galaxy is given by the col- 
lisionless Boltzmann equation, which after being integrated 
over all velocities yields the Jeans equations relating the 
potential and density of a system to its velocity dispersion 
(Binn ey" fc Tremaine 1994] l. In particular, the velocity dis- 
persion profile is important in determining the stationary 
state of the system. For equilibrium stationary states we 
found the velocity dispersion for a spherically symmetric 
stellar system, using Eq. 1201 . to be 



M 1 

a p(r a ) 



dx p{x) 



dF(x) 
dx 



(23) 



where we have assumed that = Vg and v r = ve = 0. We 
have used similar expressions to generate initial conditions 
for proto-galaxies in order to study the influence of scalar 
fields on the transfer of angular momentum during collisions 
l |Rodriguez-Meza et al 200 1) . 

The distribution function of the isotropic mod- 
els can be calculated using the Abel transform 
(Binn ey~k Tremaine 1994| . First, we define £ = —(a/M)E, 
E the total energy, and ^{r a ) = — (a/M)$jv(r a ), and we 
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obtain r as a function of \P by inverting. The mass density 
is 



p(r) 



4ty 



M x3/2 



(f) 



f(£)y/2(* - £)d£ . 



(24) 



The inversion of this equation is done using the Eddington 
formula and we obtain the distribution function 



r£ d 2 p 



The projected velocity dispersion for the isotropic sys- 
tem is given by jBinney fe Tremaine 1994} 

r dr 



4(R) = 2 



E(7?) 



p(r)vf 



R? 



(26) 



where E(i?) = 2 dr p(r)r/y/r 2 — R? is the projected den- 
sity. This also gives us the cumulative surface density 

S(R) =2tt f Z{R)RdR. (27) 
Jo 



3 POTENTIAL-DENSITY PAIRS 

In the past some potential-density pairs have been stud- 
ied within the framework of Newtonian gravity. The Jaffe 
( Jaffc 1983 ) and Hernquist density models jHernquist 1990| l 
have been generalised by Dehnen 11993^ by introducing a 
free parameter (7) that determines particular density models 
for galaxy description. Another density profile of interest, 
obtained from N-body cosmological simulations, has been 
proposed by NFW lINavarro et al 1996-9711 . In this section 
we will consider both models in order to study the influence 
of scalar fields on spherical systems. 



3.1 Dehnen's profile 

We use the family of density profiles for spherical halos and 
bulges of galaxies proposed by Dehnen (I1993II . 



p(r) 



(3-7)M 



4tt 



(28) 



ri(r + a) 4 "T ' 

where a is a scaling radius and M denotes the total mass. 
The Hernquist profile corresponds to 7 = 1 and Jaffe's to 
7 = 2. 

Solving the Poisson equation without a scalar field, i.e. 
for Newtonian gravity, the potential that corresponds to the 
density of Eq. |28J is (IDehnen 199311 



M 



1 

"2- 7 



.2-7 



In- 



; 7/2 

; 7 = 2 



(29) 



Due to the influence of the scalar field, however, ip(r) does 
not represent the total gravitational potential of the system. 
The effective Newtonian potential is now $jv = \hoo{r) 
determined by Eq. 1131 . Thus, using Eq. 12811 for the density, 
we obtain for r a < R a 



a ^ 

M* N = 



dx 



m(x) 



(30) 



r a /X a 



sinh(r a /A a ) 
+ r a /X a ^ 



and for r a > R a 

m(R a ) 



a a 



e- r °/ A ° 
r a /\ a 



-p(Ra 



(31) 



Using Eq. J2gJ one finds explicitly the functions m, p, 
and q: 



. (25 r (ro)=(3 " 7) /° dK (^)^ 



p(r a ) = (3-7)A„ 



= (3 - 7) 



dx 



(x + 1) 



4—7 



\r a + lj 
sinh(a;/Aa 



3-7 



„3-7. 



3-7 



Fi(3-7,4-7;4-7,-r a ) 



+-^0^^2^(5-7,4-7;! 



6A 2 5 - 7 



- 7, -r a ) + ■ 



(r») = (3- 7 )A a 



1-7 



dx 



-x/\ a 



(X + 1) 



4-7 



= (3 - 7) 
1 



A a 



7?r 7 2Fi(2-7,4-7;3-7,-fl a ) 



3-7 



2-7 

i?a" 7 2fi(3-7,4-7;4-7,-^a 



+ 



1 



^ 7 2^i(4 - 7, 4 - 7; 5 - 7, -R a ) - 



2Aa 4 - 7 



A 



2-7 



a 2 — f 



rr 7 2 J Fi(2-7,4- 7 ;3- 7 ,-r a ) 



„3-7. 



+ 



3-7 
1 1 

2Al~ 7' 



Fi(3-7,4- 7 ;4-7,-r a ) 



4 — 7 



Fi(4-7,4- 7 ;5-7,-r a ) - 



where 2-F1 (a, b; c, z) is a hypergeometric function. 

For r a < R a the integral in the first term of Eq. I3U1 can 
be evaluated, giving for the standard Newtonian potential 



M 



1 

2-7 



\Ra+l) Vr + 1/ 



2-7 



hi 



-Rq+1 r a +l 



7/2 

7 = 2, 



which for R a — *• 00 reduces to Eq. 12911 . 
The circular velocity is, for r a < R a , 

-r a /\ a 



a 2 



m(r a ) 



+ V 



1 + 



+ (l - y coth(r /A )) 



\ J r /A P(ra) 

Aa / T a j Aa 

sinh(r- a /A a ) 



(36) 



and for r a > R a 



m(R a ) 



M 



+ V 1 



r a /\a 



P(Ra) 



q(r a ) 



(37) 



A a / T a j A a 

The first right -hand-side term in both Eqs. I36H and 
l|37|l is the contribution from tj> and their second term 
is the contribution from the scalar field. The latter term 
depends on R a , \ a and rj = a/X a , rather than on a 
alone, which is typical for forces exerted on point systems 
l |Fischbach fe Talmadge 1 999 1 . Therefore, the stringent con- 
straints on a found for point masses at local scales are less 
severe here, since the above terms can still yield a large am- 
plitude, even when a is small. One sees that for r ^> R, \, 
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Figure 1. Rotation curves (a) and velocity dispersion curves (b) 
for the Hernquist model. The upper curves were computed with 
A a = 1 and 77 = 1, while the lower curves were obtained for A a = 1 
and 77 = 0, i.e., without a scalar field. 



the influence of scalar field decays exponentially, recovering 
the standard Newtonian result. 

The velocity dispersion can be computed by substitut- 
ing Eqs. I36H and 1371 into Eq. 12311 . and this determines 
completely the initial conditions for spherical galaxies. Fig. 1 
shows the circular velocity and the velocity dispersion curves 



(s r 



for the Hernquist model (7 = 1) for 77 = 1 (up- 



per curves) and 77 = (lower curves), i.e., the upper curves 
are computed with the scalar field taken into account and 
the lower curves without a scalar field. In the limit A — > 00 
one recovers the Newtonian result with a modified G as ex- 
pected. The computations were done with R — > 00. 



3.2 NFW profile 



JV-body cosmological simulations (Navar ro et al 1996-97fl 
suggest a universal density profile for spherical dark halos 
of the following type: 

pcrit Sc 



p(r) = 



(38) 



(r/r 3 )(l + r/r 3 y ' 

where r s is a scale radius, S c is a dimensionless constant, 
and p cr it is the critical density (for closure) of the Universe. 
This profile is the same as Dehnen's with 7 = 1 (Hernquist 
model) for r<r s , but it differs at other radii. 

This halo density profile was obtained for cosmological 
simulations in which the Newtonian theory of gravity was 
used. Because we are considering A <C H~ l , then all predic- 
tions of such cosmological simulations must be taken into 



account; this is because in the present formalism for r A 
one gets the Newtonian theory. 

Solving Eqs. @ and © with the density profile 13811 . 
we get for r a < R a , 



M 



$jv(r a ) 



-p(r a ) + 



sinh(r a /A a ) 
r a /X a 



q(r a ) 



and for r a > Ra 



m(R a ) e 



(39) 



(40) 



where the reference length a — r s . The circular velocity is 
given by Eqs. I|36|l and l|37^ where the functions m(r a ), o(r a ), 
p{r a ), and q(r a ) are given by 



m(r a ) = 
o{r a ) = ■ 

p(r a ) = 



l + r a 
1 



+ ln(l + r„)-l 
1 



1 + r a 1 + R a ' 



(41) 
(42) 



l + r a 
+ cosh(l/A c 

— sinh(l/Aa 



h(r„/A a ) 



Chi,i±^)-Chi(i. 



(43) 



<l( r a) 



A a 



,-ra/Aa 

l + r a 



-Xa 



+e 



l/Xa 



Ei 



1 + Ra 

1 + r, 



A a 



-Ei - 



1 + Ra 
Xa 



(44) 



where Ei(z), Shi(«), and Chi(z) are the exponential, hyper- 
bolic sine and hyperbolic cosine integrals, respectively. 

In Fig. 2 we plot the circular velocity and its dispersion 
for the NFW profile. The values of A a and 77 are the same as 
in Fig. 1. The upper curves correspond to the case in which 
we take into account the scalar field, while the lower curves 
were computed for standard Newtonian gravity. We took the 
size of the distribution of mass to be R — 20 r s . 

The influence of the scalar field is to enhance, especially 
in the inner regions, both the velocity curves and dispersion. 
These effects are a little more pronounced in the Hernquist 
model than in the NFW model. 



4 DISCUSSION AND CONCLUSIONS 

We have found potential-density pairs of spherical galactic 
systems within the context of linearised scalar-tensor theo- 
ries of gravitation. The influence of massive scalar fields is 
given by <j>, determined by Eq. 115H . General expressions 
have been given for circular velocities and dispersions of 
stars in the spherical system; by stars we mean probe parti- 
cles that follow the dark halo potential. Specifically, these re- 
sults were used to find potential-density pairs for the generic 
Dehnen density parametrisation, as well as for the NFW pro- 
file. In general the contribution due to massive scalar fields 
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Figure 2. Rotation curves (a) and velocity dispersion curves (b) 
for the NFW density profile. A a = 1 and upper (lower) curve was 
computed with rj = 1 (rj = 0). 



is non-trivial, see for instance Eqs. I36137L and interest- 
ingly, forces on circular orbits of stars depend on the am- 
plitude terms R a , X a and r\ = a/\ a . This means that even 
when local experiments force a to be a very small num- 
ber ( |Fischbach fc Ta lmadge 1999), the amplitude of forces 
exerted on stars is not necessarily very small and may con- 
tribute significantly to the dynamics of stars. Alternatively, 
one may interpret the local Newtonian constant as given by 
+ instead of being given by {<t>)~ 1 (1 in our con- 

vention). In this case, the local measurement constraints are 
automatically satisfied, and at scales larger than A one sees 
a reduction of 1/(1 + a) in the Newtonian constant. If this 
were the case, then the upper curves in figures 1 and 2 would 
have to be multiplied by 1/(1 + a). 



In the past, different authors have used point solutions, 
Eqs. © and © , to solve the missing mass problem encoun- 
tered in the rotational velocities of spirals and in galaxy clus- 
ter dynamics ISanders T984llEclKardt 19931 . These models 
were used as an alternative to avoid dark matter. Indeed, 
for single galaxies one can adjust the parameters (a, X) 
to solve these problems without the need for dark matter. 
However, these models do not provide a good description 
of the systematics of galaxy rotational curves because they 
predict the scale A to be independent of the galactic lu- 
minosity, and this conflicts with observations for different 
galactic sizes l |Aguirre et al 200H unless one assumes vari- 
ous A's, and hence various fundamental masses, ra, one for 
each galaxy size. Such a particle spectrum is not expected 
from theoretical arguments; it represents a considerable fine 



tuning of masses. This criticism would also apply to our 
models. The purpose of the present investigation was, how- 
ever, not to present a model alternative to dark matter in 
order to solve the missing mass problem, but to compute the 
influence of scalar-field dark matter distributed in the form 
of a dark halo. This contribution, together with bulges and 
disks, will give rise to a flat velocity curve, as in Newtonian 
mechanics. Yet, the influence of the scalar field from STT 
yields some profile modification, especially in inner regions, 
as is shown in figs. 1 and 2. 
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